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Abstract 

A new method to extract the density of partition function zeroes (a continuous function) from their distribution 
for finite lattices (a discrete data set) is presented. This allows direct determination of the order and strength 
of phase transitions numerically. Furthermore, it enables efficient distinguishing between first and second order 
transitions, elucidates crossover between them and illuminates the origins of finite-size scaling. The efficacy of the 
technique is demonstrated by its application to a number of models in the case of Fisher zeroes and to the XY 
model in the case of Lee- Yang zeroes. 
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1. Introduction 

Phase transitions are phenomena common to 
a wide range of disciplines, from physics to biol- 
ogy, economics and even sociology. Examples in- 
clude neural networks, protein folding, melting, 
magnetism, stock market crashes and the decon- 
fincment transition in the early universe. 

In statistical physics, in particular, one is inter- 
ested in the determination of the location, order 
and strength of phase transitions. First order tran- 
sitions involve the coexistence of two or more dis- 
tinct phases and are characterised by a discontinu- 
ity in the first derivative of the free energy corre- 
sponding to the finite amount of energy needed to 
transform from one phase to the other. For temper- 
ature driven transitions, this discontinuity is the 
latent heat Ae and is a measure of the transition 



strength. For second order transitions there is no 
such coexistence. Instead, thermodynamic quanti- 
ties such as the correlation length and specific heat 
diverge. Such divergences are characterised by crit- 
ical exponents v and a, which is a measure of the 
strength of the transition. 

The computational approach to the study of 
phase transitions consists of two steps - the gath- 
ering of data in the form of a Monte Carlo simu- 
lation followed by a numerical analysis of appro- 
priate quantities. The first of these is restricted to 
systems of finite size. Phase transitions, however, 
require an infinite number of available states for 
their occurrence. The second step in the numerical 
approach is an extrapolation to infinite volume. 
Traditional techniques involve the finite-size scal- 
ing (FSS) study of thermodynamic functions. An 
increasingly popular alternative approach is, how- 
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ever, the use of zeroes of the partition function. 

Let t = T/T c — 1 be the reduced temperature 
and h the external field. For a d-dimensional sys- 
tem of linear extent L, the FSS of the j th complex 
partition function zero (for large j) is given by [1] 



t,(L)^(j/L d ) 



1/vd 



or 



hj(L) ~ (j/L d ) 



{d+2- v )/2d 



(1) 



(2) 



Here rj is the anomalous dimension and ft, = in 
the first formula (where the zeroes are called Fisher 
zeros) while t = in the second (where the zeroes 
are Lee- Yang zeroes). The standard approach to 
FSS of zeroes is to fix the index to j = 1 and ex- 
tract an estimate for the critical exponents from 
a range of lattice sizes. It has, however, long been 
known that using more than one index could pro- 
vide more information. However, since (1) and (2) 
are inexact, this has been prohibitive. In particu- 
lar, the extraction of the density of zeroes (a con- 
tinuous function) from their (discrete) distribution 
for a finite and numerically accessible lattice has 
been considered prohibitively difficult or even im- 
possible [2]. In recent years, however, there have 
been some attempts to overcome the difficulties in- 
volved [3] . In view of the increasing importance at- 
tached to this approach, we suggest an appropriate 
way this should be done [4] . 



2. Density of Zeroes 



The partition function for finite L is Zl(z) oc 
Y[j (z — Zj(L)), where z is an appropriate function 
of temperature or field. We assume the zeroes, Zj, 
are on a line impacting on to the real axis at the 
critical point, z c . Parameterising zeroes on this line 



by Zj 



+ rj cxp (tip) we may define the den- 



sity of zeroes as <7z,(r) — L d J2j 3( r ~ r j(L))- The 
cumulative distribution function of zeroes is then 
Gh{r) = ^ gh{s)ds which is j / L d if r E (r,-,r,+i). 
At a zero we assume the cumulative density is given 
by the average [5,6] 



In the thermodynamic limit and for a phase tran- 
sition of first order this integrated density of zeroes 
is, in fact [5], 



Goo(r) = .9oo(0)r 



(4) 



so that the density is non- vanishing at the real axis. 
The slope at the origin in (4) is related to the latent 
heat in the Fisher case (or magnetization in the 
Lee- Yang case) via [5] 



>(0) oc Ae 



(5) 



For a second order transition the corresponding 
expressions for Fisher and Lee- Yang zeroes are [7] 



Goc(r) oc r 



2-a 



or Goo(r) oc r 2d/(d+2-r,) 5 (6) 



G L {r 3 ) = (2j-l)/2L d 



(3) 



respectively. 

Thus while the scaling behaviour of the position 
of the first few zeroes in the complex temperature 
plane can be used to identify v via (1), the density 
of zeroes gives the strength of the transition. A plot 
of Gh(rj) against Tj(V) should (%) go through the 
origin, {it) display L- and j- collapse and (Hi) re- 
veal the order and strength of the phase transition 
by its slope near the origin. 

In (3), rj may be taken to be the imaginary 
part of the position of the j th zero. Equating 
(3) to (6) in the second order Fisher case, gives 
r-{L) - L- 1 /". This is the usual FSS formula 
for fixed index Fisher zeroes. Similarly, in the 
Lee- Yang case, one recovers the fixed index FSS 
formula hj (L) - L"(< i + 2 -'7)/ 2 . Also, equating (3) 
to (4) gives rj(L) ~ L~ d , explaining the usual 
identification of v with 1/d for a first order tem- 
perature driven phase transition. Therefore, tra- 
ditional FSS emerges quite naturally from this 
density approach. 



3. Applications 

To demonstrate our approach, we perform fits 
to the cumulative density of zeroes for a number of 
different models in statistical physics and in lattice 
field theory [4] . Here we present the results for the 
Fisher zeroes for two models, one from each field. 
Furthermore, we discuss the Lee- Yang zeroes of the 
two-dimensional XY model. 



2 



Allowing for first or second order behaviour, the 
cumulative density should behave as 



G(r) = a ir a2 + a 3 



(7) 



where we also allow for an additional parameter 03 
which should be zero for a good fit. In fact, a non- 
zero value of 03 indicates the absence of a phase 
transition, for, a% > means the zeroes have al- 
ready crossed the real axis (the situation in the bro- 
ken phase) while 03 < means the zeroes have not 
yet reached the real axis (the symmetric phase). 
For Fisher zeroes, a first order phase transition is 
indicated if ai ~ 1 for small r, in which case the 
latent heat is proportional to the slope a\. A value 
of (ii larger than 1 signals a second order transi- 
tion whose strength is given by a = 2 — 02- Note 
that a can be measured directly using this method 
while traditional FSS only allows the measurement 
of the ratio a/ v. 

The d = 2, q= 10 Potts Model: The first six Fisher 
zeroes for the two-dimensional 10-state Potts 
model for lattice sizes L — 4-64 are listed in [8] . A 
traditional FSS analysis applied to the first zero 
for large lattices provides evidence for v = 1/d 
and hence a first order phase transition. However, 
the determination of the lattice size above which 
FSS sets in is, by necessity, somewhat arbitrary. 
Indeed, when one extends the analysis to higher 
index zeroes one finds that when corrections are 
ignored, no two-parameter fit gives an acceptable 
result. 

Our analysis of the density of zeroes begins with 
Fig. 1. The excellent data collapse for various L 
and j indicates that (3) is the correct form for the 
density of zeroes. Fitting (7) to the L = 16-64, 
j = 1-4 data points gives 02 = 1.10(1) and a 3 = 
0.00004(1), indicating a first order phase transi- 
tion. Fixing a 3 = 0, a 2 — 1 and applying a single- 
parameter fit close to the origin yields a slope cor- 
responding to latent heat Ae = 0.698(2) which 
compares well with the exact value 0.6961. 

The d = 4, Abelian Surface Gauge Model: This is 
a model dual to the d = 4 Ising model, which, up 
to logarithmic corrections has mean field critical 
exponents [9]. One therefore expects the surface 



CD 



0.03 



0.02 



0.01 




0.04 



Fig. 1. Distribution of zeroes for the d = 2, q = 10 Potts 
model (for L = 16-64) which has a first order phase tran- 
sition. The symbols X, +, A,o, □, and O, correspond to 
j = 1, 2, 3, 4, 5, and 6, respectively. 

gauge model also to be characterised by mean field 
exponents with a = and v = 1/2. 

The first two Fisher zeroes for lattices of size 
L = 3 to 12 are listed in [10] where a conventional 
analysis applied to the first index zero yields the 
best estimate of v = 0.469(17) from the two largest 
lattices. Inclusion of the smaller lattices worsens 
the fit driving v away from 1/2. Also, a bimodal 
structure in the energy histograms appears as a 
spurious indication of a first order transition [10]. 

A fit of the data to (7) yields a-i incompatible 
with unity (see Fig. 2), with a fit near the origin 
yielding 02 = 1.90(9). This corresponds to a = 
0.10(9), compatible with zero. Note that only the 
region near the origin in Fig. 2 is of interest. The 
slope there is not compatible with a first order 
transition and the bimodal structure of the energy 
histograms observed in [10] can only be an unex- 
plained finite-size effect. 

The d = 2 XY Model: Here we demonstrate that 
the density technique is also applicable in the Lee- 
Yang case. The first few Lee- Yang zeroes for the 
d = 2 XY model with L = 32-256 were deter- 
mined for various temperatures in Ref. [11]. Fig- 
ure 3 depicts the distribution of these zeroes for 
(3 = 1/ksT = 1.113 which is the critical point 
(here ks is the Boltzmann constant). From (6), and 
with 77 = 1/4 (the value of the anomalous dimen- 
sion in this model), one expects Gi(rj) ~ r* 6 ^ 15 ', 
where rj = hj. A three-parameter fit to (7) gives 
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Fig. 2. Distribution of zeroes for the four-dimensional 
Abclian surface gauge model which has a second order 
transition. The symbols □ and O correspond to the j = 1 
and j = 2 index zeroes, respectively. 
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Fig. 3. Distribution of Lee- Yang zeroes for the 
two-dimensional XY model at (3 = 1.113 (for L = 32—256). 
The symbols X , +, A, o, and □ correspond to j = 1, 2, 3, 4, 
and 5, respectively. 

= indicating that criticality has indeed been 
reached at this temperature. A two-parameter fit 
now yields a 2 = 1.06(1), compatible with expec- 
tation. One notes that logarithmic corrections are 
present in this model as shown in [11,12]. 



4. Conclusions 

We have presented a new method to extract 
the (continuous) density of zeroes from (discrete) 
finite-size data and demonstrated how this can be 
used to distinguish between phase transitions of 
first and second order as well as to measure their 
strengths. The method meets with a high degree of 



success in statistical physics and lattice field theory 
and lends new insights into the origins of finite-size 
scaling. 
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